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Abstract 

Diffusion of information, behavioral patterns or innovations follows diverse pathways depending 
on a number of conditions, including the structure of the underlying social network, the sensitivity 
to peer pressure and the influence of media. Here we study analytically and by simulations a general 
model that incorporates threshold mechanism capturing sensitivity to peer pressure, the effect of 
‘immune’ nodes who never adopt, and a perpetual flow of external information. While any constant, 
non-zero rate of dynamically-introduced spontaneous adopters leads to global spreading, the kinetics 
by which the asymptotic state is approached shows rich behavior. In particular we find that, as a 
function of the immune node density, there is a transition from fast to slow spreading governed by 
entirely different mechanisms. This transition happens below the percolation threshold of network 
fragmentation, and has its origin in the competition between cascading behavior induced by adopters 
and blocking due to immune nodes. This change is accompanied by a percolation transition of the 
induced clusters. 

There are remarkable analogies between the social contagion of information, behavioral patterns or 
innovation and some physical or epidemic spreading processes, where global phenomena emerge through 
the diffusion of microscopic states W- All evolve in networks with nodes characterized by relevant state 
variables, and links that represent direct interactions between nodes. In biological systems epidemics 
are driven by binary interactions that lead to the emergence of simple contagion phenomena [1]. Social 
diffusion processes are usually characterized by complex contagion mechanisms, where node states are 
determined by comparing individual thresholds with all neighbor states HISHH]- This property, capturing 
the effect of peer pressure and commonly assumed in social spreading phenomena 13 HQ], has consequences 
on the dynamics and the final outcome of the social contagion process. Moreover, the theoretical approach 
to these systems has much in common [niiiii], which greatly helps us to understand their behavior. 

Models employing threshold mechanisms mostly focus on cascading phenomena where, under some 
circumstances, a macroscopic fraction of nodes in the network is converted rapidly due to microscopic 
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Figure 1. (color online), (a) Numerical simulation of a general threshold model over an empirical 
network, with adoption threshold (p = 0.2, rate of spontaneous adopters p = 0.0005 and fraction of 
blocked nodes r = 0.1. The network is an ego sample of Facebook friendships with size N = 96 and 
average degree z = 10.63 [IS]. Susceptible nodes adopt spontaneously with rate p or after a fraction 
(p of their neighbors has adopted, while blocked nodes never adopt, (b) Schematic illustration of the 
spreading process. At t = to node P spontaneously becomes an adopter, ‘infecting’ nodes a and b. When 
Q adopts, it induces the adoption of nodes c — f. Nodes inside the ellipse constitute an induced cluster 
of adoption. 


perturbations. This approach is motivated by earlier social theories mm and has been implemented by 
Watts in an elegant model of cascading behavior |S|. Watts showed that a global cascade (occupying a 
macroscopic fraction of the network and induced by local perturbations) can occur due to the interplay 
between network structure and individual thresholds. He further identified the phase with a non-zero 
probability of global cascades in the space {p, z) of the average threshold p of nodes and the average 
degree z of the network. 

While the relevance of this model is indisputable |61[II1[I1H23|, its limitations become clear from real 
social spreading data. The Watts model focuses on the (instantaneous) emergence of global cascades 
triggered by single local perturbations, while there are empirical examples where threshold mechanisms 
do play a role yet global adoption phenomena emerge through other scenarios. In reality global adoption 
is often not induced by microscopic perturbations but by a larger fraction of people [22] . Moreover, deci¬ 
sions of individuals depend on external impulses arriving from mass media or advertizing |24j . resulting 
in a perpetual stochastic perturbation. In addition, there are individuals entirely reluctant to adopt. 
Furthermore, the Watts criterion for macroscopic adoption is purely deterministic, coded in the network 
structure, threshold distribution and perturbation site - it does not concern time, which is clearly a 
feature of empirical stochastic processes of adoption spreading. 

Here we present a general threshold-driven model of social contagion phenomena that captures various 
spreading scenarios, ranging from cascading behavior to dynamically evolving non-explosive patterns, and 
sheds light to the different kinetics behind them (Fig. [^. Motivated by empirical observations [2S|, we 
extend Watts’ threshold model by considering blocked nodes immune to social influence and discuss their 
effect on cascade formation. In addition, we introduce spontaneous adopters with a constant rate, and 
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Figure 2. (color online). Frequency Fg of global cascades as a function of node threshold (p and network 
average degree z at an intermediate time t = 100, for varying p and r. (a-b) As p increases, the region 
that allows global cascades of adoption grows in size. Its boundary is well approximated by a cut-off in 
the fraction of adopters p as calculated by Eq. ([^ . (c-d) Conversely, an increasing fraction of blocked 
nodes shrinks the global cascade regime. Dots show the boundary of this regime according to Eq. Q. 
Simulations correspond to an ER network with N = 10^ and are averaged over 10"^ realisations. 


present approximate analytical and numerical results regarding our model. In particular, we study how 
the kinetics of spreading changes for an increasing density of blocked nodes. We aim at the simplest 
possible but sufficiently general extension of earlier threshold models [SI with a minimal set of 

states and transitions necessary to describe various real scenarios of social spreading phenomena. The 
introduction of further states, secondary adoption, or other decision-making mechanisms is left as a 
further challenge, since our aim here is to model generic cascades of primary adoption. 

In Watts’ threshold model |5], all nodes are initially in a susceptible state 0, except for a single 
adopter seed in state 1. The process evolves as each node with degree k changes its state from 0 to 1 if a 
fraction p of its neighbors have adopted before. Since nodes cannot change their state after exposure, the 
system evolves towards a state where no further adoptions are possible. The emergence of a global cascade 
depends on the degree distribution p^ of the network, the distribution p^ of individual thresholds, and the 
initial seed. The condition for a global cascade is the existence of a percolating component of vulnerable 
nodes with thresholds Q < p <l/k (who need one adopting neighbor before exposure) connected to the 
seed. This percolating vulnerable tree is quickly converted after adoption of the seed and may trigger 
further adoption of stable nodes with thresholds p > Ijk (who need more than one adopting neighbor to 
adopt). Assuming an Erdos-Renyi (ER) random network [26] and a single adopter seed, there is a phase 
boundary in {p, z)-space encompassing a regime where global cascades occur (Fig.j^a). The properties of 
this cascading regime have been investigated for the case of heterogeneous thresholds, different network 
topologies Enm, and variable seed size mm- 

Empirical studies, however, support the intuition that some individuals in society may refuse to adopt 
technological innovations for various reasons - due to another favorite product, aversion towards a firm. 
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or some criticism on principle |25j . Such individuals will never be exposed, irrespective of the state of 
their neighbors m- To consider this behavioral pattern in our model, we block the adoption of a fraction 
r of randomly selected nodes in the network. These nodes do count when their neighbors consider the 
decision to adopt, and thus will make it harder for neighbors to fulfil the threshold criterion. 

Accordingly, the original Watts model corresponds to r = 0, while for r > 0 the phase diagram 
changes. Even in the presence of blocked nodes, macroscopic spreading is still determined by the static 
criterion of the existence of a global vulnerable cluster, and thus a generating function technique [6] can 
be applied [^[28]. Assuming a single threshold (j) and an ER network with average degree z, the condition 
for the emergence of a macroscopic cascade is. 


(1 - r)e 




— z = 0, 


( 1 ) 


with kc = [!/</*] • to the factor 1 — r, the introduction of blocked nodes shrinks the region in 

(0, z)-space where global cascades develop, in good agreement with numerical simulations (Fig. [^c and 
d). 

While blocked nodes hinder the spreading process, there are reasons other than social influence that 
could motivate individuals to adopt a social pattern, like external influence from mass media. This 
spontaneous adoption has been studied theoretically by introducing a given density of adopters at the 
outset of the Watts model [22] . However, spontaneous adopters may get active at any time during a real 
social contagion. Thus we include a stochastic dynamics where a susceptible node may become adopter 
with rate p at any time, irrespective of the status of its neighbors. 

Considering both extensions, we have a threshold-driven dynamics with three node states: blocked, 
susceptible and adopter (FigQ. At the outset, all nodes are susceptible except for a fraction r that 
remains blocked. At each time step of the simulation, a randomly selected, susceptible node i adopts 
spontaneously with probability p, otherwise it adopts if at least a fraction (j) of its neighbors has already 
adopted. If r = 0 and p > 0 all nodes will eventually adopt (Figj^a), following a kinetics reminiscent of 
the approach to a unique ground state in a physics system. On the other hand, if we introduce quenched 
randomness and stochastic perturbations {r,p > 0), our model allows various temporal regimes and a 
transition from rapid to slow spreading dynamics. 

Our threshold model can be studied analytically by extending the framework of approximate master 
equations (AMEs) for monotone binary-state dynamics developed by Gleeson [T5l - H7l |29|, where the 
transition rate between susceptible and adoption states only depends on the number m of neighbors that 
have already adopted. We ignore topological correlations by considering a configuration-model network 
with degree distribution p^. and average degree z. We describe a node by the property vector k = (fc, c), 
where fc = 0,1,... is its degree and c = 0,1 its type, i.e. c = 0 is the type of the fraction r of blocked 
nodes, while c = 1 is the type of all nodes that may adopt with threshold (j). Moreover, pk is a joint 
distribution giving the probability that a randomly selected node has property vector k. Assuming 
independence between degrees and types, pk = rpk for c = 0 and pk = (1 — r)pk for c = 1. 

The rules of our model are condensed in the probability that a k-node will adopt in a small 

time interval dt, given that m of its neighbors are already adopters, where. 




p if m < k(j) 
1 if m > k(j) 


Vm and k > 0, 


^see Supplemental Material (SM) for details on the analytical treatment of the model. 


( 2 ) 
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Figure 3. (color online). Numerical simulations and analytical approximation of the threshold model 
for 2 ; = 7 and (j) = 0.2 (continuous and dotted lines, respectively), (a) Fraction of adopters p as a function 
of time for varying p and hxed r. (b) Time evolution of the normalized adoption density p/(l — r) for 
different values of r and fixed p. (c) Final relative density /p°° as a function of r, for both spontaneous 
and induced adopters (i = 0,1, respectively), (d) Final fraction of spontaneous adopters p^ as a function 
of r. (e) Normalized maximum speed of spreading Prn,/[p(^ — ?')]i calculated from the derivative of p{t). 
Shaded areas signal the regime of slow contagion r > Vy,. Curves correspond to N = 10'^ and are averaged 
over 10^ realisations. 


with = 0 \/k,m and F’(o.i).o =P (for blocked and isolated nodes, respectively). The dynamics of 

adoption is well described by an AME for the fraction Sk,m(0 of k-nodes that are susceptible at time t 
and have m = 0,..., k adopting neighbors [13 [m HH] , 

^k,m — Ak^m^k.m r7l)sk,m “t” Ul l)5k,m—(3) 

where, 

Q _ SkPk ~ ru)Fk,m'Sk,m /.x 

To reduce the dimensionality of Eq. ([^ we focus on p{t), the fraction of adopters in the network, and 
^(t), the probability that a randomly chosen neighbor of a susceptible node is an adopter. We consider 
the ansatz Sk,m = for m < kcj) with the binomial distribution ~ 

leading to the condition b = /3s(l — i^). Then, the AME system is reduced to the pair of ordinary 
differential equations (see SM), 


p = h{v,t) - p, 


with initial conditions p(0) = i^(0) = 0. Here, 


h{v,t) = {l-r) ft + il- ft)^Pk 

k mykcj) 


(5a) 

(5b) 


( 6 ) 
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and, 


(7) 


= {1-r) /t + (1 -/i)^-Pfe Bk-i.mii') 


m>k(p 


with ft = 1 — {1 — p)e A linear stability analysis of the reduced AME system recovers the cascade 
condition for p = 0 0 (see SM). Moreover, the fraction of adopters p{t) obtained by solving Eq. 0 is 
in considerable agreement with numerical simulations (Eig. [^a and b). Since susceptible nodes adopt 
spontaneously with rate p, their fraction po{t) in the network is approximated by, 


Po=P [ (1 - r - p)dt. (8) 

Jo 


where p{t) follows Eq. 0 (Fig.|c and d). We denote its counterpart pi = p — po ss the fraction of 
induced adoptions, i.e. vulnerable and stable adopters. 


For p > 0 the dynamics has a trivial asymptotic state with a final fraction of adopters p°° = 1 — r, 
however, the kinetics of the model depends on the parameters. We first focus on the frequency Fg 
of global cascades (i.e. adoption reaching at least 20% of susceptible nodes [S]) and its behavior in 
{(p, 2 ;)-space for varying p and r. For fixed t and p > 0, there is a region where global cascades occur 
(Fig. I b) that can be compared with the asymptotic cascade regime found for p = 0. The boundary 
of this regime is well approximated by Eq. Q for p = 0 and by Eq. ([^ for p > 0. By continuously 
introducing spontaneous adopters the global cascade regime expands, meaning that macroscopic adoption 
is eventually possible for systems with any degree and threshold. Even in the absence of a percolating 
vulnerable component in the network, a growing number of spontaneous adopters induces local cascades 
that merge due to triggered stable adoptions and finally form a giant component. This behavior is 
consistent with empirical observations in the online spreading of communication technologies |25j . 

The kinetics of spreading may change by introducing many blocked nodes. As r (and thus random 
quenching) increases the adoption process slows down (Fig. [^b). In this dynamics nodes change state 
in two ways: i) via spontaneous adoption (a slow process for small p), or ii) via induced adoption 
by fulfilling the threshold condition, which may lead to fast cascading behavior. For small r induced 
adoptions dominate spreading (Fig. |^c) and p grows rapidly towards p°°. On the other hand, for large 
r adoption slows down since stable nodes have more blocked neighbors and it is difficult to fulfil their 
threshold condition. This slow regime is mostly driven by spontaneous adoption, as evidenced by the 
relatively large asymptotic fraction of spontaneous adopters pjf (Fig. id). 

Taking the ER network as example, a giant component of susceptible nodes can only exist for r < 
r* = 1 — 1/z [Sni El] . Then, a relevant question is whether regimes of fast and slow spreading are 
separated by a characteristic value rx < r*. One possibility is to define rx as the value that maximizes 
p^ and for which ~ with the final fraction of induced adopters. For z = 7, ^ = 0.2 and 
p = 0.0005 we have « 0.7 and r* = 0.857, meaning that slow spreading occurs even in susceptible 
networks that are not fragmented. The slow regime is further characterized by the lowest possible value 
in the maximum spreading speed, pm ~ p(l ~ ?'); corresponding to the rate of spontaneous adoption at 
the beginning of the dynamics (Fig. [^e). In other words, the time series p{t) has an inflection point for 
r < Tx and is concave for r > rx. 


To better understand the kinetics of the crossover between spreading regimes around rx, we finally 
focus on the size distribution P{s) of induced clusters, i.e. connected components of adopters disregarding 
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Figure 4. (color online), (a) Size distribution P(s) of induced clusters in the regimes of fast {r < rx) 
and slow (r > rx) spreading, at several stages in the adoption process. In early times, P{s) is unimodal 
and qualitatively similar in both regimes. As t increases, the distribution becomes bimodal only for 
r < Tx, indicating the presence of global cascades, (b) Asymptotic size distribution Poo(s) of induced 
clusters, after t = 5000 and for varying r. For r < Vx global cascades may still develop and make Poo(s) 
bimodal. As r increases, the distribution becomes unimodal and global cascades disappear. Simulations 
correspond to z = 7, (j) = 0.2, p = 0.0005, N = 10"^, and are averaged over 10^ realisations. 

spontaneous adopters (Fig. 0 b). For early times P(s) includes small induced clusters only, indicating that 
a larger fraction of spontaneous adopters is crucial for global spreading in the absence of a percolating 
vulnerable component (Fig. |^a). However, for late times the behavior of P(s) differs between regimes: 
in the regime of rapid spreading the distribution becomes bimodal due to the appearance of a global 
cluster of induced adopters, while in the slow regime (r > rx) it remains unimodal until the end of the 
dynamics. Overall, the crossover between regimes seen globally in the speed of spreading (Fig. is 
accompanied by an underlying, percolation-type transition revealed by the asymptotic size distribution 
Poo(s) (Fig.|4]b). Indeed, in the asymptotic limit t —>■ oo and as r increases, this distribution stops being 
bimodal at r « 0.74 ~ rx for the studied ER case. 

The peculiarity of this dynamic percolation transition of induced clusters is that, in contrast to static 
percolation problems, it is not known a priori which node will participate in the process, as any unblocked 
node may become a spontaneous innovator. By analyzing the properties of this transition, we find a 
critical percolation point at rc ~ 0.738 with the set of exponents (3 = I.l, 7 = 1.0, r = 2.5, and iz = 3.1, 
which are rather close to the mean field values El 

Our aim in this paper has been to provide a general dynamic model of social spreading phenomena 
that accounts for various kinetics. Our model is designed such that it: (a) is driven by threshold 
mechanisms capturing the role of social pressure, and (b) concerns temporal aspects of the emergence of 
global cascades. We generalized Watts’ threshold model [H] with mechanisms of spontaneous adoption 
and complete reluctance to adoption, in order to further understand the temporal behavior of spreading 

^For notation and the related theory see D. Stauffer and A. Aharony, Introduction to Percolation Theory (Taylor and 
Francis, London, 1994); except that the exponent —{1/u) is defined here not by the scaling with the linear dimension of 
the system but with its size. 
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phenomena. We have shown that, outside of the cascading regime of the Watts model, there is possibility 
of global contagion mediated by spontaneous adopters. However, the speed of spreading depends strongly 
on the density of blocked or immune nodes. For a small fraction r of blocked nodes, few spontaneous 
adopters enable the formation of large clusters by initiating cascades. For large r, spreading slows down 
as it is dominated by spontaneous adopters and only small cascades are generated. Our intrinsically 
dynamic model is able to describe various scenarios of real social contagion as well as the crossover 
between them, and shows a novel percolation transition of induced clusters. This model has not only the 
potential to explain observational data |25j but, with appropriate fitting, may help identify the character 
of spreading processes at an early stage, hinting in this way at possible measures to improve adoption 
performance. Moreover, it is possible that the consideration of blocked nodes will help understand a 
diversity of spreading phenomena, including related seismic or neural processes. 

RZ acknowledges support from FP7 MULTIPLEX Grant No. 317532, GI from the Academy of 
Finland, and JK from H2020 FETPROACT-GSS CIMPLEX Grant No. 641191. 
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Supplemental Information for 

Kinetics of social contagion 

Z. Ruan, G. Ifiiguez, M. Karsai, J. Kertesz 


In this Supplemental Information we present two different schemes of analytical treatment for the 
dynamical threshold model introduced in the main text. First, in Section]^ we solve the model exactly 
for the case of a single innovator seed (p = 0), but allowing for a non-zero fraction of blocked nodes to 
be present in the system (r > 0 ). We do this by extending the generating function approach provided 
by Watts [5]. Second, in Section]^ we provide an approximate solution for the general case of multiple 
innovator seeds introduced periodically in time (p > 0) and a non-zero fraction of blocked nodes (r > 0). 
This solution is based on the approximate master equation (AME) formalism introduced by Gleeson 

[iMiziisn]- 


A Generating function approach for p = 0 

Here we consider the case p = 0, r > 0, i.e. there are no spontaneous adopters in the network (except 
for a single seed), but a fraction r of blocked nodes hinders the spreading process. Similarly to the 
Watts case [6l [22], the phase diagram can be explored as a static, percolation problem. Our goal is 
the condition for the existence of a giant vulnerable component 0, which can be found by using the 
generating function approach [55] • Suppose pk is the probability that a randomly chosen node has k 
connections, and pk is the probability that a node with degree k satisfies the condition 1/k > (j), where (/> 
is the threshold of the node. Since a fraction r of nodes is blocked, the probability of a randomly chosen 
node being unblocked and satisfying the threshold condition is /Ofe(l ~ f)- These nodes are vulnerable, 
since they may adopt only if at least one of their neighbours has already adopted. Thus the probability 
of node with degree k being vulnerable is PkPk{^ — f). The corresponding generating function is, 

Goix) = '^PkPki^ - r)x’". (9) 

k 

Another generating function we are interested in is Gpx), which generates the degree distribution 
of a vulnerable node b that is a random neighbour of a given node. The probability of choosing node b 
with degree k is kpk/z, thus we have, 

G,(^) = ^ ggM, ,10) 

z z 

where z = J2k is the average degree of the network. 


II 




To calculate the size distribution of the connected components consisting of vulnerable nodes, we 
consider two other generating functions, 


Hoix) = 

71 

(11) 

Hl{x) = y^^WnX", 

(12) 


n 


where Qn is the probability that a randomly chosen node belongs to a vulnerable cluster of size n, and 
Wn is the probability that a random neighbour of a given node belongs to a vulnerable cluster of size n. 


First let us calculate Hi{x). Notice that a random graph below percolation can be regarded as a 
tree-like structure (since the probability of containing a loop scales as N~^ and is negligible for large N 
[28|1. Under this assumption, a random neighbour of a given node a can be in one of several states: it 
may not be vulnerable; it may be a vulnerable node with no edges (other than the one connecting to 
o); it may have one edge connected to another component, two edges connected to two components, etc. 
Then Hi{x) takes the form, 


Hi{x) = WnX” = U(x") 

n 

= Pi+xGi{Hi{x)) 

= \-Gi{l) + xGi{Hi{x)), (13) 


where Pi is the probability that a random neighbour of a is not vulnerable, meaning that the node 


belongs to a vulnerable component with size n = 0. According to Eq. (10), 01(1) is the probability that 
a random neighbour of a is vulnerable, thus we have 1 — Gi(l) = Pi. The term xGi{Hi(x)) comes from 
the fact that Hi{x) satisfies a self-consistency condition [55]. Furthermore, Ho{x) can be calculated in a 
similar way, 


Ho{x) = 1 - Go(l) -f xGo{Hi{x)). 


(14) 


The average vulnerable cluster size is (n) = Hq{1). Using Eqs. (10), (13) and (14) we have, 


(n) — Go(l) -I- 


(Go'(l))^ 

z-G"(l)- 


(15) 


This expression is similar to the one found for the Watts model |6]. However, we have to use Eq. [^for 
Gq. Thus (n) diverges for. 


Gq (1 ) =y^ k{k- 1)(1 - r)pkPk = z. (16) 

k 

When Gq( 1) < z all vulnerable clusters are small since (n) is finite, while for Gq( 1) > z there is a giant 
vulnerable cluster percolating throughout the system. 


In principle we may solve Eq. (16) for arbitrary pk and (j) (or pk). We now restrict ourselves to the 
special case of a Poisson distribution pk (corresponding to an Erdos-Renyi random network) and constant 
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Figure 5. Cascade window for different values of r. Dashed lines separate two different phases: in the 
inner region global cascades may happen, but not outside. 
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(j) (i.e. all nodes have the same threshold). In this case, pk = e ^z^/kl and pk satishes. 


Pk = 


1 , k < kc 
0 , k > kc 


where kc = [!/</>) • Substituting the expressions oi pk and pk into Eq. (16) we get. 


k^2 


(17) 


(18) 


an expression that can be solved numerically. Notice that there are 3 parameters here, z, r and (j). We 
may, for example, vary any two of them (z and (j)) and solve for the third one (r). 


As it is shown in Figj^ the cascade-allowing phase in {(j), z) space depends on the value of r. Each 
dashed line encloses the region in which the cascade condition is satished for a given r. As r increases, the 
cascade window shrinks in both (j)- and 2 :-axes, as blocked nodes hinder the formation of a giant vulnerable 
component. The lower boundary of each phase diagram is mainly constrained by the connectivity of the 
network (since most nodes satisfy the threshold condition because of their low degree) and can be 
estimated by a mean-field approximation. In average each node has z connections, and among these z 
neighbors there are zr blocked nodes, meaning that the effective degree of a node is z — zr. The condition 
for the existence of a giant vulnerable cluster is an average effective degree larger than 1. Thus we have 
z = 1/(1 — r), determining the lower boundary of the cascade window. For r = 0.5 we have z = 2, which 
agrees well with calculations shown in Figj^ 


In Figj^ (a) we show the critical r value coming from Eq. ( [T^ as a function of the average degree z 
and threshold (p (only positive values are viable solutions). According to Eq. (16), the calculated critical 
r values determine a contour surface below which there is a giant vulnerable component; otherwise all 
vulnerable clusters are small. Figj^ (b) shows the projections of this surface on the (</, z), (</, r) and 
(z, r) planes. Notice that in the projection on the (z, r) plane, r as a function of z first increases and then 
drops down. The reason is that when the average degree z is very small or very large, the corresponding 
largest vulnerable cluster is small, i.e., only a few blocked nodes are needed to destroy it. 
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Figure 6. (a) Critical r value allowing for the emergence of giant adoption clusters as a function of (p 
and z. (b) Projections of subplot (a) on the (^, z), (^, r) and (z, r) planes. 

B AME formalism for the general case 

B.l Stochastic binary-state dynamics 

Here we extend a general rate equation formalism for stochastic binary-state dynamics as developed 
recently by Gleeson [ISHnlES]. In a stochastic binary-state dynamics, each node in the network can 
take one of two possible states (susceptible or adopter in the language of innovation adoption) at any 
point in time, and switching state randomly with probabilities that only depend on the current state 
of the updating agent and on the states of its neighbors. This general definition includes our threshold 
model as a special case. Such formalism considers configuration-model networks, that is, an ensemble 
of networks specified by the degree distribution pk but otherwise maximally random (i.e. pairs of stubs 
are connected uniformly at random, so that in the limit iV —>■ oo of very large network size there are no 
degree-degree correlations or clustering). 

For each node to describe all of their relevant properties we introduce a vector k = (fc,c), where 
fc = 0 , 1 ,... fcjvf is the degree of the node and c = 0 , 1 a dummy variable that labels its ‘type’, i.e. any 
other property that characterizes the node apart from its degree. In the case of our threshold model, 
c = 0 is the type of the fraction r of blocked (or immune) nodes, while c = 1 is the type of all non- 
blocked nodes that adopt with threshold p. The maximum value /cm is tuned to approximate the degree 
distribution with any level of accuracy. Any pair of nodes with identical values of k are considered 
equivalent in this level of description, forming a node class with the same average dynamics. Moreover, 
Pk can be generalized to the joint distribution giving the probability that a randomly selected node 
has property vector k (i.e. degree k and type c). If blocked nodes are chosen randomly among all nodes, 
like in our model, then pk = rpk for c = 0 and pk = (1 — r)pk for c = 1 . 
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In the language of innovation adoption, the dynamics of a node is determined by the number m = 
0 ,1 ,... fc of its neighbors that have already adopted when the node is deciding whether to adopt or not. 
During a small time interval dt, a susceptible node (in state 0) adopts with probability Fk^mdt, while an 
adopter (in state 1) becomes susceptible with probability R^^ rndt. The functions and Rk.m, known 
as infection and recovery rates, respectively, determine the temporal evolution of the node class k. In the 
particular case of threshold models, a so-called monotone dynamics, i?k,m = 0 Vk, m (since no adopters 
become susceptible again). As for Fk,mj the rules of spontaneous and threshold adoption imply. 


I P if m < kd) 

, Vto and k > 0, (19) 

1 if m > k(j) 

that is, a node adopts the innovation either spontaneously with rate p, or with probability 1 if its number 
of adopting neighbors equals or exceeds the integer threshold $ = [ kef]. Blocked nodes (c = 0) have an 
infection rate of = 0 Vfc, m, while for isolated nodes {k = 0) F(^Q i'f Q = p. In other words, blocked 

nodes never adopt, and isolated nodes can only adopt spontaneously. 


Let us now turn to the rate equations for our threshold model, called AMEs in the formalism by 
Gleeson. We denote by Sk,m(i) the fraction of k-class nodes that are susceptible at time t and have m 
adopting neighbours. Therefore, the fraction of agents with property vector k that are adopters at time 
t is Pk(t) = 1 — J2m=o 'Sk,m(t)) and the fraction of adopters in the system is p{t) = X]k^'kPk(t)- Here, the 
sum over classes means a sum over all degrees and types, i.e. Ec Assuming a monotone 

dynamics {Rk^m = 0), the AMEs for Sk,m can be written as [TSIITBII^ . 


dsk. 


dt 


= -Tk,mSk,m - /3s(fc - rn)sk,m + Psik -m + l)sk,m-l, 


( 20 ) 


where m = 0,..., fc, Sk,-i = 0, Fk,m follows Eq. (19), and /3s(t) (the rate at which edges between pairs 
of susceptible nodes transform to edges between a susceptible agent and an adopter) is given by. 


Psit) = 


SkPkZm(^ ~ w)Ek,mgk,m(t) 

EkPkEm(^-™)'Sk,m(t) 


( 21 ) 


If at time t = 0 there is an infinitesimally small seed for the adoption process (i.e. p(0) = 0), the initial 
conditions for Eq. (20l are Sk.m(O) = Bk,m{P), with Bk^m a binomial factor, 

BkM=(^\p^{l-pf-^. ( 22 ) 


The solution Sk,m{t) of the AME system in Eq. (201 provides a very accurate description of the 
dynamics of our model, yet its dimension is (kM + 1)(^m + 2). Therefore, the number of equations 
to solve grows quadratically with the maximum degree. Fortunately, the AMEs for our model can be 
mapped to a reduced-dimension system with a derivation similar to the one used by Gleeson in the case 
of the Watts threshold model [51 HI] . 


B.2 Reduced-dimension AMEs 


To reduce the dimension of Eq. (201, we need to consider system-wide quantities that are more ag¬ 
gregated than Sk,m- One of them is the probability that a randomly chosen node is an adopter, p[t) = 
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1 —Sk,m(i)j i-6- fraction of adopters in the network. The other one is the probability that a 
randomly chosen neighbour of a susceptible node is an adopter, v{t) = 


We start by proposing an exact solution for the AME system in terms of the following ansatz, 

Sk,m(t) = Bk,m[v{t)\e~^^ for m < k<j) and c = 1, (23) 


and = Bk^m{v) for c = 0, where Bk^m follows Eq. (22|. The meaning of the ansatz in Eq. (23) 

is quite intuitive and considers two processes. First, a susceptible agent with degree k and m adopting 
neighbours is connected to m adopters with the binomially distributed probability Second, 

for m < kcj) a susceptible node does not fulfill the threshold rule and can only adopt spontaneously 
with probability since the system is progressively been filled by adopters. Considering these two 


processes as independent we end up with the product in Eq. (23). Finally, since blocked nodes do not 


adopt and are distributed randomly over the network, S(fe,o),m is determined only by a binomial factor. 


The next step is to insert the ansatz (23) into the AME system (20) and derive a set of differential 


left-hand side of Eq. (20)) we get. 


equations for the aggregated quantities p and v. Taking the time derivative Sk,m of Eq. (23) (i.e. the 

(24) 


Sk.m 


m k — m 

V \ — V 


V-p\ 


Then, we use the threshold rule (19) for to < k4>, the ansatz (231 and the binomial identity, 

Bk^ra—liy^ — , p B}^ Yniy') ^ 

V K — TO -I- 1 


in the right-hand side of Eq. (20) to obtain 


--fhc.mSk.m - f^s(k - TO)sk,m + f^sik - TO -|- l)sk.m-l = 


\ — V 

—p + j3s [ m — k -\ -TO 


Sk,7 


Equating Eqs. (24) and (26) as in the AME system (20) leads to 


i> = Psil - v), 


(25) 


(26) 


(27) 


a condition on v so that the ansatz (23) is a solution of Eq.(20). This differential equation has the 


initial condition i^(0) = p(0) = 0, obtained by evaluating Eq. (23) at t = 0 and comparing with the 


expression i?fe,m(0), which corresponds to an infinitesimally small amount of initial adopters randomly 
distributed among k classes. Furthermore, by assuming a (yet to be determined) function g{v,t) such 


that i> = g{v,t) — u, Eq. (27) reduces to 


Ps = 


g{v,t) - V 

\-v ■ 


(28) 


Now, we consider the following general result derived by Gleeson in [T^ (Eqs. (F6)-(F10) therein), 

^Pi^^{k - m)sk,rn = z{l - vf, (29) 
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with z = ^Pk the average degree in the network. Eq. (291 is valid for functions Sk,m and v that satisfy 


Eqs. (20) and (27) respectively, for any Ek,m and random initial conditions on Sk.m and v, and is thus 


applicable in our case. Our goal here is to use Eq. (29) to find an expression for g{v) and therefore write 


the differential equation (27) explicitly. Noting that the left-hand side of Eq. (|29[) is the denominator in 


the definition (21) of /3s and that = 0 (i.e. blocked nodes do not adopt), Eq. (21) gives 


Ps = 


1 — r 

1 

z(l-j.)2 


ik-m)s^k,i),ra + J2 Pk ^ ^ ^)^(A;,l),r; 

k m<.k4> k m>k(p 

- m)sk,m - - w)s(fc, 0 ) 


-{l-r){l-p)'^Pk ^ {k-m)s(k,i),r. 

k m<k(f) 


(30) 


where we have written pk explicitly as pk = rpk for c = 0 and pk = (1 — r)pk for c = 1 , in order to 


stress the dependence on r. Then, we insert the ansatz (23) (with its special case S(k,o),m = Bk^m{v) for 
blocked nodes), as well as the identities (k — ra)Bk^rn(v) = k{l — v)Bk-i^rn{k') and J2m<k<p Bk-i,m{k') = 
1 - Em>fc 0 Bk-i,m{v) to obtain, 


Ps = 


1 - V 


( 1 -r) 


l-(l-p)e ^’*-f(l-p)e ^*'^-Pk X! Bk-i^miy) 


m^kip 


(31) 


A comparison of Eqs. (28) and (31) gives us the following expression for g{v,t), 


g{’^,t) = {l-r)ift + {l-ft)'^-Pk 


(32) 


m>k(p 


where we define /t as /t = 1 — (1 — p)e Thus, the AME system (20) gets reduced to the differential 


equation i> = g{v,t) — v, with g{v,t) given explicitly by Eq. (32). 


Even though the equation v = 5 ( 1 ^, t) — is closed and in this sense equivalent to Eq. (20), we can 
also derive the corresponding equation for p, since we are mainly interested in the temporal evolution of 


the fraction of adopters in the network. From the definition of p and Eq. (20) we have, 

P ~ ^ ^ ^ ^ '5k,m ~ ^ ^ P\^ ^ ^ Ak y,.jSk,m 

km km 

-I- /3s ^ Pk ^ [(A: - m)sk,m - {k - m + l)sk,m-i] , (33) 


where the second term in the right-hand side telescopes to zero. Then, we use an algebraic manipulation 
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similar to that of Eqs. (30) and (31) to obtain, 


,m 5k, m — (1 I P ^ ^ Pk ^ ^ ,m E Pk X! '*('=.!) ,ri 

fc m<.k(f) k m'>k(p 


k m 


= (l-r) I l-(l-r)(l-p)^pfc S(fc,i),m ) - P 

k m<.k(l) 


= (1 - 0 I /t + (1 - ft)YP>= E ^k,m{v) j - p. 

k m^kcj) 


(34) 


In this way, Eqs. (33) and (34) can be rewritten as p = h{v,t) — p, where. 


h{v,t) = (1 - r) I /t + (1 - ft)YPk E I • 

k m'>k<p 


(35) 


Joining all of these results, the AME system (20) gets reduced to the system of two ordinary differ¬ 
ential equations. 


p = h{v,t) - p, 
i' = t) - i/. 


(36a) 

(36b) 


with the quantities g{v,t) and h(i',t) given explicitly by Eqs. (32) and (35). 


The system (36) can be solved numerically to obtain p{t) and thus characterise the temporal evolution 
of the adoption process. Let us further separate the fraction of adopters as p{t) = po{t) + pi{t), where 
po and Pi are the fractions of innovators and induced adopters, respectively. Now consider the identity 




(37) 


where Ps{t) is the fraction of non-blocked, susceptible nodes that can eventually adopt, either spon¬ 
taneously or not. Since such suceptible nodes adopt spontaneously at a rate p, the rate equation for 


innovators is po = PPs- Then, with Eq. (37) we obtain. 


Po{t) = P [ [^ - r - p{t)]dt, 
Jo 


(38) 


which can be calculated explicitly with the numerical solution of Eq. (36). 


B.3 Cascade condition for p = 0 


We may also use the AME formalism in the case p = 0 (no spontaneous adoption) to derive a cascade 


condition similar to Eq. (16), as has been done previously for the Watts model [5S]. First we observe that 
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by assuming p = 0, /t = 0 and thus we remove the explicit time dependence in Eq. (36), i.e. g = g{v) 
and h = h{v). Then the sum over m in g{iy) can be rewritten as, 


k k 

= (1 - r) ^^ Bk-i,m{i') = {i-r)'^-pk'^Bk-i,m{’^)fik,m), (39) 

k m>k4> k m—0 


where. 


f{k,m) = 


0 if m < k(j) 
1 if m > k(j) 


(40) 


and /(0,0) = 0. Eq. (40) is the so-called response function of the monotone dynamics in our model, i.e. 
a function that activates when a node with degree k and m adopting neighbours fulfils the threshold 
condition and therefore may adopt. 


Let us now perform a linear stability analysis of the reduced AME system (36) around the equilibrium 


point (p*,j^*) = (0,0), corresponding to a total lack of adoption. If (p*,i^*) is unstable, then any small 
perturbation (like a single node adopting at t = 0) can drive the system out of equilibrium and create a 
global cascade of adoption where p > 0, that is, a system where a non-vanishing fraction of nodes has 


adopted in the limit N ^ oo. Since the equation u = g{i/) — n is closed, the stability of Eq. (36) is 
determined by the stability of this equation at = 0. Then, according to linear stability theory, a local 
instability exists at = 0 if. 


dv 


bM- 


> 0 , 


(41) 


that is, if d^g{0) > 1, where d^, denotes derivative with respect to v. 


We can write d,^g explicitly by inserting the definition of the binomial factor into Eq. (32) 

dg 


k-1 


du 


= -Pk fik, m) 


m—0 


—l,n 

du 


with 


dBi. 


k—l.m 


fc-i\ 

V 

m J 


\i-vf 


— {k — 1 — m)^'"(l — v) 


fe —2—ml 


(42) 


(43) 


dv \ m 

Then we analyse terms in the sum over m at the equilibrium = 0. For m = 0, dj/i?fc-i,o(0) = 1 — fc, 


but f{k,0) = 0 for 0 > 0 and thus the term in Eq. (42) is zero. For to = 1 we have d,^Bk-ip{0) = k — 1. 


Finally, for to > 1 we get dvBk-i^m{^) = 0. Overall, the instability condition (41) gets reduced to. 


k 

This condition defines the area in (0, z)-space where global cascades may develop, for a given value of r. 


By comparing Eqs. ( [T7| ) and (40) we see that f(k, 1) = p^, and so the boundary described by Eq. (16) is 
recovered by Eq.([44|. In other words, the cascade regime for p = 0 is described accurately by both the 
generating function approach and the AME formalism. 
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